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“TELE 9 OTssroucs - COMPLES ** INTEGER*4 power routine 
File OTSPOWCJ.MAR Edit: 5BL 1004 
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COPYRIGHT (c) 1978, 1980, 1982, 1984 
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FACILITY: Language support library - user callable 
"ABSTRACT: 


COMPLEX base to INTEGER*4 power. 
Floating overflow can occur. 

Undefined eppgnent tet ten can occur if 
base = (0.,0.) and exp <=0 


VERSION: 0 
aereet? 
Jonathan M. Taylor, 29-jun-77: Version 0 
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Edit history for version 0 
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HISTORY ; Detailed Current Eds t History MTHRTL.SRCJOTSPOWCJ.MAR; 1 


-SBTTL HISTORY ; Detailed Current Edit History 


anes SbeRtOR © NTH SSIGNAL JMT 5-0CT-77 
han Ms rome e &° JACKET to MTH FLAG, JACKET. TNH 17-July-78 
- - Fix ug givin ng aiv de fault, or wrong results for 

ve He | _ clean up comments. 


; 1-001 = Change at Von number tg 1 god Hr. UNDE XP 
; to ATHSK_UNDEXP. 58S 7-D 
Add te” the PSECT directive. S 22-DEC-78 
Declare externals. SBL 1 
= Use general mode addressing. "SBL 30-Nov-1981 


; eistery fer Version 0 0 OTSsPowc4 
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te INTEGER*4 power ro shen 16-SEP-1984 01:56:54 VAX/VMS M v04-00 
NS ¢ : gre eer1o8e 91338300 REJOTSP 


$6 ~SBTTL DECLARATIONS 


65 : INCLUDE FILES: 


o8 : EXTERNAL SYMBOLS: 
0 e 


IGNAL 3; Math error routine 
- ; COMPLEX division routine 


> MACROS: 


: EQUATED SYMBOLS: 
base = 


: OWN STORAGE: 


PSECT DECLARATIONS: 


-PSECT _OTSSCODE PIC,SHR,LONG,EXE ,NOWRT 


3000000009000 09 09 09 09 09 09 09 SININI INI NINN NINO AOAAAAO 
AUIE WN SO OONAULSWN ("OO OONOUS WN NOue 
. 
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4 ; base input formal - by-value 
exp = 12 3; exponent intpu formal - by-value 


; program section for OTS$ code 
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«SBTTL OTSSPOWCJ = COMPLEX*8 ** INTEGER®4 


p++ 
; FUNCTIONAL DESCRIPTION: 


COMPLEX result = COMPLEX base ** INTEGER*4 exponent 
The COMPLEX result is given by: 


base exponent result 
any >0 PRODUCT (base * 2**i) where 
is each non-zero bit in 

exponent. 

(0., 0.) <=0 Undefined exponentiation. 

not (0., 0.) <0 PRODUCT (base * 2**i) where 
i is each non-zero bit in 
texponent:. 

not (0., 0.) =0 (1.0, 0.0) 


tt- 
Floating overflow can occur. 
Undefined exponentiation occurs if base is 0 and 
exponent is 0 or negative. 
CALLING SEQUENCE: 
Power.wfc.v = OTSSPOWCJ (base.rfc.v, exponent.rl.v) 
INPUT PARAMETERS: 
The base input parameter is standard FORTRAN COMPLEX. 
The exponent input parameter is a signed longword integer. 
Both input parameters are CALL BY VALUE. 


IMPLICIT INPUTS: 
NONE 

OUTPUT PARAMETERS: 
NONE 

IMPLICIT OUTPUTS: 
NONE 


FUNCTION VALUE: 
COMPLEX base ** INTEGER*4 exponent 
SIDE EFFECTS: 
SIGNALS SS$_ARITH with floating overflow hardware code if 
floating overflow. 
SIGNALS MTHS_UNDEXP (82 = * UNDEFINED EXPONENTATION') if 
base is 0 and exponent is 0 or negative. 
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UNDEF INED: 
ASHQ 


—— GOOD 000000000 0000000900 09 60 09 09 09 09 09 09 SI NIN NINN NN NIAAA AAAOAOOUIUIUIN 
ROOD NA NE WIN OOD NAME WN (OOO NAUE WIN $C OONAULS WN $0 OONOAOUE WO OOnNO 
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MTHRTL.SR €10 OTSPOWCJ.MAR; 1 
-ENTRY OTSSPOWCJ, “M<R2,R3,R4 »R5,R6> 


g eigapte integer jy ten 
MOVQ base(AP), R2 ; R2/R35 gets COMPLEX bas 
MOVL exp(AP), "RG ; R4 = Longword exponent. 
BGEQ 1$ ; R4 = | exponent : 
MNEGL R4, RG 
BBCC #0, R4, EVEN ; branch if even and clear low bit 
MOVa ; RO/R1 = initial result 
ROTL #-1, RG, RG 3; R4 = unsigned_exponent / 2 
BEQL Done 3; done if exponent was 1 


N 
BRB SQUAR1 else use rest of exponent 


MOVD 1 0 ; RO/R1 = initial result 
ROTL i, R4, R4 3 R4 = unsigned. exponent / 2 
BNEQ SQUA : branch if exponent not 
TSTF 2 3 exponent was 0, text RP(base) 
BNEQ DONE 3; done if non-0, answer is 1. 
TSTF ; IP(base) better not be zero 
BNEQ DONE 3; it isn't return 1.0 

#15, #1, RO ; return RO = reserved operand 
MOVZBL aMTHsK GNDEXP -(SP) 3 FORTRAN —e number 
CALLS , G*ATHS$SIGNAL 3 convert to 32-bit condition code 

; and SIGNAL MTH$_UNDEXP 

RET 
ASHL #-1, R4, R4 ; R& = treduced exponent: / 2 


= square current base 


MULF3 R2, R3, RS : R5 = tmp = RP(base)*IP (base) 
MULF R2, ¢ ; R2 = RP(base)**2 
MULF R35, R : R35 = IP(bas oeee 
SUBF R3, R2 $ Ré = RP(base) #* - IP(base)**2 
ADDF3 = R5. RS, R3 : = 2*RP(base)*IP(base) 
LBC R4, SQUAR ; branch if next exponent bit is 0 
= partial result * current power of base 
MULF3 Re: R3, RS : RS = tmp = RP(part) * IP(base) 
MULF 3 RO = RPtpart t) * RP(base) 
MULF3 RI, ~ R6 3; R6 = IP(part) * IP(base) 
SUBF e 3; RO= RPC Dart) #RP base) -1P (pa rt) *1P(base) 
MULF R2, R1 ; R1 = 1P(part) * RP(base) 
ADDF R5, R1 3; Ris SP yer cae cecee sie \geceseertbees! 
ASHL = #1, RG, R4 : R4 = {reduced expon / 2 
BNEQ SQUAR1 ; loop if more cnpanent «hits left 
TSTL dy yl ; test exponent sign 
BGEQ POWCJ : done if positive 
TSTF 3; test RP(result) 
BNEQ RECIP : if non-0, OK to take reciprocal 
TSTF R1 > RP(result) was 0, test IP(result) ;tt- 


7 
| 
erage SrGDGRES cAMPERG*stoveeGaNeL™ CREE AUSSEG) WAMAME SSB IOS anes" cy 


AF «(13 13 eecip BEQL UNDEF INED 3; undefined (0.040.011) ** <n 
7E 30 7D 15 mMOVQ RO, -(SP) 3 Second ers. getr is divisor 
, Bee 70 16 MOVD S*#1, -(SP) 3 push on stack 
00000000 ' GF 04 FB 0088 17 CALLS #4, 6*OTS$DIVC H Bort s reciprocal 
008F 18 POWCJ: jtt- 
04 O08F 19 RET 3; result in RO/R1 
0090 $9 
0090 1 - END 
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The_working set Limit was 900 pages. 
3043 bytes (6 pages) of virtual memory were used to buffer the intermediate code. 
There were 10 pages of symbol table space allocated to hold 13 non-local and 1 local symbols. 
- source lines were read in Pass 1, producing 11 object records in Pass 2. 

pages of virtual memory were used to define 0 macros. 


BASE = 00000 0% | 
DONE 5 R 4 
EVEN 1F R 1 
EXP = oc 
MTHS$SIGNAL teeerere x 0 
MTHSK_UNDEXP eerereee x 0 
OTSSDTVC eereeere = § x 0 
OTSS$POWCJ 00000000 RG 01 
CJ QOOOOOBF R 3 

RECIP 0000 ore R 1 
SQUAR 00000041 R 01 
SQUAR1 Bopoeege R 4 
UNDEF INED 0000031 R 1 

Pewee ooececocecs 4 

! Psect synopsis ! 

pow mewn meme wom + 
PSECT name Allocation PSECT No. Attributes 
s Oe s 00000000 0.) 00 ¢ OO.) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
_OTSSCODE 00000090 ¢«¢ 144.) O01 ¢ 1.) PIC USR CON’ REL LCL SHR EXE RD NOWRT NOVEC LONG 

tem nem em mmo answerer escecmee 
; Performance indicators ! 

Phase Page faults CPU Time Elapsed Time 
Initialization 3 00:00:00.09 00:00:00.89 
Command processing 12 i SS 00:00:06.31 
Pass 7 00:00: 3-69 00:00:04.26 
Symbol table sort 0 BF 380: 88- 1 00:00:00.10 
Pass 2 54 00:00:00.46 00:00:01.69 , 
Symbol table output 2 4s Sh B88) 20 «Bs 
Psect synopsis output 2 00:00:00. 00:00:00.0 
Cross-reference output 0 4 Se ge 
Assembler run totals 289 0:00:01.76 00:00:15.35 


Geowcceesoeweeeeeocoeccccce$ 


Geoweeeeeoeeecsooceesoeeser$ 


Macro Library name Macros defined 


-$255$DUA28: CSYSLIBISTARLET.MLB; 2 0 
0 GETS were required to define 0 macros. 


There were no errors, warnings or information messages. 


TS 


OTs CJ 
VAX-11 Macro Run Statistics 
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MACRO/ENABLE=SUPPRESSION/D1I SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:OTSPOWCJ/OBJ=OBJ$:OTSPOWCJ MSRC$:OTSPOWCJ/UPDATE=(ENH$:OTSPOWCJ) 
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